Nodal Lines in the Cranked HFB Overlap Kernels 
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Norm overlap kernels of the cranked Hartree-Fock-Bogoliubov states are studied in the context 
of angular momentum projection. In particular, the geometrical distribution of nodal lines, i.e., 
one dimensional structures where the overlap kernels possess null value, is investigated in the three 
dimensional space defined by the Euler angles. It is important to know the distribution of these 
nodal lines when one attempts to determine the phase of norm overlap kernels. 



The mean field approximation (MFA) is successful in 
the description of nuclear systems 0, @, , to which 
the independent particle motion picture can be applied 
0- 0- Deformation and superfluidity (superconduc- 
tivity) in nuclei are good examples |3j that are explained 
by MFA. Although some essential correlations are taken 
into account through MFA in an efficient manner, it fails 
to incorporate other important features such as symme- 
tries and relevant higher order correlations, which may be 
necessary for studying nuclear structures further Q, y| . 

Within MFA, physical quantities such as energy and 
quadrupole moment are calculated as expectation values, 
but there is room for further quantisations. For exam- 
ple, angular momentum is not conserved when the ro- 
tational symmetry is spontaneously broken by deformed 
mean fields. As a consequence, the corresponding many- 
body state has the form of a wave packet @. In the 
cranking model , which describes a rotating mean field 
in a semi-classical manner, the mean field state is written 
as 
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(The rotational frequency of the deformed mean field is 
denoted by to, while the quantum numbers for total an- 
gular momentum and its magnetic quantum number are 
expressed as I and M, respectively [22J-) IV'mfa(^)} usu- 
ally has large fluctuations in its probability distribution 
IC^I 2 . The width of the fluctuations becomes larger at 
higher spin [6j], which means |^mfa(w)) possesses a more 
averaged character. As a consequence, a simple applica- 
tion of the cranking model faces difficulty, for instance, 
in the description of the mixture of two different states 
in band crossing regions |8| • 

A quantum mechanical description of collective rota- 
tion is given in the generator coordinate method (GCM) 
with a choice of the Euler angles (O = (a, 0, 7)) 



as generator coordinates. This method corresponds to 
angular momentum projection (AMP) |Tcj . Rotational 
symmetry spontaneously broken by the deformed mean 
field can be restored by superposing the states that point 
various orientations \ip(ti)) = -R(0)|t/>mfa) ■ ( R(Q) is a 
rotational operator in the three dimensional space. See 
Eq.@ below.) A state with symmetry restoration |<E>) is 
schematically expressed as 



$) = / dn/(n)wn)>, 
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with /(ft) being the weight function that is determined 
by the variational principle. This state |<fr), a superpo- 
sition of the infinite number of degenerate states, corre- 
sponds to a Nambu-Goldstone mode [ljj in a finite sys- 
tem. 

In our previous works, we have numerically performed 
AMP in attempts to describe high-spin states, such as 
tilted rotation for high-if states in 178 W 0] and the 
wobbling motion in the multi-band crossing region in 
182 Os [l3j. Through these works, we found it difficult 
to perform numerical calculations of AMP in some situ- 
ations. Such situations are observed to occur, for exam- 
ple, in the calculations of AMP applied to cranked HFB 
states at high spin (angular momentum ranging from 10 
to 20h) . A typical problem is negative values of the prob- 
ability \C ! M 1 2 for / being odd integer, although the prob- 
abilities for even / seem reasonably calculated |23f ■ (See 
FigE) 

We first suspected the coarse mesh employed in the 
numerical integration causes this difficulty, because the 
integration in the AMP formula (See Eq.Q below) is 
discretised in numerical calculations with respect to the 
Euler angles. The appropriate mesh size can be esti- 
mated by using the uncertainty principle AiAft ~ fi. 
The constraint value ( J\) = J and the fluctuation AI are 
roughly related as AI ~ J/2 at high spin [y, [Tj| . Then, 
Afl = 2h/J. For example, Aft ~ 6° for J = 20h. The 
appropriate number of mesh points are therefore about 
60 for a, 7, and 30 for (3. Even with this mesh size, how- 
ever, the problem was observed to occur in performing 
AMP when high spin states are considered (i.e., J > \QTi). 
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FIG. 1: Probability distributions of angular momentum using two types of prescription to determine the sign of the norm 
overlap kernels. The solid line represents the probability obtained by the old method while the dashed line is by the new 
method. The constraint value on angular momentum is J = 14h. The left (right) panel shows even-/ (odd-/) components. 
Note the different scales in the first and second panels. 



On the other hand, when the angular momentum is low 
(J < 10h), AMP is successfully performed, as expected. 

We then suspect that the poor determination of the 
sign of the norm overlap kernels, which will be given in 
Eq.(J7J, may cause this problem. These kernels are nec- 
essary in the AMP procedure, in which the state in the 
projected space is expressed as 
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The coefficients g T K arc determined by a variational equa- 
tion 



E ( H mk - E'Nl^) g l K = 0. 



(4) 
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Here, P\ IK denotes the AMP operator, and matrices 
N[ 1K and H T MK are defined as, 
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A useful relation to note here is |C^| 2 = Nj^ M . The 
measure in the three-dimensional Euler space dSi is de- 
fined as dasin/3d/3e?7, and the integration intervals are 
[0, 2tt] for a and 7, and [0, tt] for (3. The rotational oper- 
ator R(Si) is given as 

R(Si) = cxp(—iaJ z ) exp(—i/3J y ) exp(—i~fj z ). (6) 

The quantities (V'mfaI-R(^)IV'mfa) and 
(iPmfa\H R(Sl) IV'mfa) in Eq.JSJ) are called norm and 
energy overlap kernels, respectively. They are calculated 
through formulae derived by Onishi et al [Tit ITr3 |. By 



utilising the Onishi formulae, one can calculate the norm 
overlap kernel as 



N(Si) = (V>mfa|£(0)|V>mfa) = <7(fi)VdctP(fi), (7) 

where a (Si) has values ±1. This ambiguity in sign comes 
from the square root operation. We will come back to 
this point shortly. The matrix P(Sl) is written 



p(si) = u*d\sl)u + v J 'D T (n)v, 



(8) 



where U and V are matrices defined in the general Bo- 
goliubov transformation [3( between canonical (c M c>) 
and quasi-particle annihilation and creation operators 

(A 
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The presence of a (Si) implies that N(Si) is a two- 
valued function of Si although it should be physically 
well-defined and single-valued. Therefore, it is necessary 
to choose an appropriate si gn for given Si. Two meth- 
ods have been proposed |r7lll8| to determine the sign of 
cr(Si), but due to its numerical feasibility the method of 
Hara, Hayashi and Ring is widely used. It makes use 
of the continuity of the norm overlap kernel as a function 
of the Euler angles Si. Our calculations are based on this 
method, and the procedure is explained below. 

First of all, the origin in the Euler space is assumed to 
have some certain phase. (In our case, <r(0) = +1.) Then, 
by analytical continuation, the well-defined region is ex- 
tended until the whole Euler space is covered. To achieve 
this, the Euler space is divided into domains defined as 
follows: Writing N(Si) = r(Si)e l6 ( n \ the argument 6(Sl) 
and the norm r(Si) are supposed to be continuous with 
respect to Si. Then, a point SI = Sl(a, (3, 7) belongs to a 
domain D n (n = 0, ±1, ±2, • ■ •) when (n - ±)tt < 9 (Si) < 
(n + jj)tt. Inside each domain, a (Si) is determined to be 
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+1 for even |n| and —1 for odd |n|. The boundary be- 
tween D n and D n+ i is determined by tracking the con- 
tinuity of 8(H) through the derivative information. In 
this way, N(il) becomes well-defined as a single-valued 
function in the entire domain D = [J n D n (i.e., the whole 
Euler space) . 

As mentioned above, derivatives of the norm overlap 
kernels are used to check the continuity of 9(H). A set 
of formula derived by Onishi et al. is again useful to 
calculate these derivatives |16). which are now given in 
the form of logarithmic derivatives, 



d_ 



ln(^MFA|^(^)|^MFA) 
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2 detP(O) Oil 



detP(il). 



Detailed forms for the formulae are given in Ref. 16] . 
(Sometimes the formulae are called generalised Wick's 
theorem Q.) Because the argument of the norm overlap 
kernel is written as, 



0(n)=lm Ih^mfaIPOWmfa) 



(11) 



the derivatives of 8(H) can be computed through Ea. ((TU|l . 
Note that the derivative of the argument, 8' = j^, is 
uniquely evaluated , i.e., being free from the sign ambi- 
guity, because these formulae essentially make use of the 
square of the norm overlap kernels (N(il) 2 ) to eliminate 
the ambiguity coming from a(il). 

To specify the domain D n , we consider the follow- 
ing two quantities that can be calculated from a pair 
of neighbouring mesh points displaced by small distance 
Ail, that is, (i) the difference between the two points, 
corresponding to the approximate derivative of the ar- 
gument 8(H) obtained by the two-point formula, which 
is 



8(n + An) - e(ii) 

a{n) = Ail 



(12) 



and (ii) the average of 9' obtained through Eq. l|10|) be- 
tween the two points, which is 



m ^ 8'(ii + Aii) + 9'(ii) 



(13) 



Note that only a(il) suffers from the sign ambiguity, but 
b(il) does not. We can now determine the relative sign of 
the kernel between the points il and il + Ail by checking 
the magnitude of the residue 88(0) = a(il) — 6(0). When 
the value of the residue is small, or behaves as 89(d) ~ 

- Ail 2 with il < il < il + Ail, the neighbouring 

point il + Ail is judged to be in the same domain D n 
as the reference point il. On the contrary, the point 
il + Ail belongs to the different domain if the following 
two conditions are satisfied: (i) 88 returns a significant 

deviation from the behaviour of — 9 Ail 2 ; and (ii) 
the deviation is removed after the sign a(il + Ail) is 
inverted and the resultant residue follows the behaviour 

< 3 >(Q) 



According to the results in our calculations, the above 
naive approach works for J < 10h, but it seems to fail 
when high spin states are considered (J > 10?i). A typ- 
ical symptom of the failure is seen as violation of the 
positive definite nature in probability, as illustrated in 
the right portion of Fig.l. In Fig.l, angular momentum 
components in the cranked HFB states at J = 14H are 
plotted, that is, 



W 1 = Ty(N ! ) = 



i 

E 

M=-I 



\c i m ((j 1 ) = iaK)\ 



(14) 



of 



'-Ail 2 . 



where the matrix N 1 is defined in Eq.JSJ) and the trace 
is taken with respect to the magnetic quantum number 
M. With the above method to determine the sign (solid 
lines in the figure), AMP seems to work properly when 
the even-/ distribution of W 1 is concerned (i.e., the left 
portion in the figure) , while in the right portion the posi- 
tive definite nature of W 1 is violated slightly to an order 
of magnitude of 10 -3 . One may say that such small nu- 
merical errors can be neglected when proper physics is 
extracted from the results. In fact, as far as the ground 
state rotational bands are concerned (which consist of 
only even- J components), it may be justified to ignore 
these small errors. However, we would like to investigate 
the cause of the errors in this paper because high-spin 
physics involving odd-/ states (e.g., high- if bands) be- 
comes more and more important today. 

The above problem (i.e., negative probability) implies 
that the resultant norm overlap kernels do not satisfy 
proper features as a physical quantity. A possible reason 
for this defect is the assumption of the continuity of 8(il). 
It is, in fact, singular on a nodal line, i.e., a set of points 
where the norm overlap kernel becomes zero (r(il) = 
0). The analytical continuation approach assumes the 
continuity of norm overlap kernels (more precisely, of the 
argument 8(H)), so that the method fails to work in the 
neighbourhood of nodal lines. 

To deal with the nodal lines in the context of the sign 
determination, there are basically two approaches. One 
is to find out all the nodal lines in the Euler space (from 
the information of r(il) — 0). Once the locations of the 
nodal lines are known, one can determine a(il) without 
any ambiguity because we can avoid the calculations of 
the residue 88 between two points sandwiching the nodal 
lines. The other is to calculate the residues for all the 
possible directions (that is, Ail =±Aa, ±A/3, and ±A7) 
no matter where the nodal lines are. In this case, opti- 
misation regarding to the values of the residues is made 
by fixing a(il + Ail) in a manner of trial-and-error. In 
this paper, the second approach is chosen for the sake of 
simplicity in numerical calculations. 

Although there is no need to know explicit information 
about the geometry and topology of the nodal lines in 
the present approach, we realise that such information 
are very useful to clarify the nature and source of the 
problem. Therefore, in this paper, we try to pin down 
the locations of nodal lines, and succeed to find them for 
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the first time. (Figs |2] and E{J. 

Let us first discuss the nodal lines we found, before 
the detailed explanation is given on the new method to 
determine the sign a. To obtain norm overlap kernels, 
we first perform principal axis cranked HFB calculations 
with a constraint on angular momentum (Ji) = J (as 
well as constraints on particle numbers). The pairing- 
plus-QQ force is employed with the standard choices for 
the force parameters and model space |l9( . This HFB cal- 
culation is performed in a fully self-consistent manner, so 




(A) Aa = A 7 = 10°, and A/3 = 2°. 




(B)Aq = A7 = 5°, and A/3 = 1°. 



FIG. 2: Nodal lines at J = 12h with low and high resolutions 
for the mesh size. 



that the gamma (triaxial) deformation is self-consistently 
handled. The details of the method used for this calcu- 
lation are explained in Ref.|23|. Then, as already men- 
tioned above, the overlap kernels are obtained by means 
of the Onishi formulae such as Eqs. Q and (JTDJ. The 
modified method for the phase determination, which will 
be explained later, is employed to calculate the norm 
overlap kernels. 170 Dy is chosen in this study, as a well- 
deformed rare-earth nucleus. (}3 — 0.295, A p = 0.867 




(A) Aa = A7 = 10°, and A/3 = 2°. 



1 > 




(B) Aa = A 7 = 5°, and A/3 = 1°. 



FIG. 3: Nodal lines at J = 14h with low and high resolutions 
for the mesh size. The detailed structure of (B) is shown in 
FigH 
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(MeV), A„ = 0.652 (MeV) at J = P.) The yrast 
line of this nucleus does not show any crossings at high 
spin (< 20K.) in our calculation. (The more complicated 
structure of nodal lines may appear in the cranked HFB 
states in band-crossing regions. Such states are planned 
to be studied in future.) 

Figs|21 and |21 present the nodal lines found in cranked 
HFB states at J = \2h and 147i, respectively. In the 
figures, cubes are used to show and cover the regions 
where nodal lines exist, for graphical convenience. The 
size of the cubes is the same as the mesh size used in 
AMP. It is given by AQ = L n /N n , where L a = L 1 = 2n 
and Lp = tt. The number of the mesh points are N a — 
N 7 = 36 (72) and N p = 90 (180) for the low (high) 
resolution calculations. 

The positions of nodal lines are determined in the fol- 
lowing approximation: on each face of the cube, after 
the phase determination, the norm kernels are linearly 
interpolated with respect to the Euler angles, that is, 

Minoar(fi) - C\OL + C 2 /3 + C37 + C 4 . (15) 

The coefficients Cj are determined by the exact values 
of N(Jl) at the four vertices on the face. The position 
of the nodal line on the face is obtained by solving the 
simultaneous equations of 

Re[AWar] = Im[iVlincar] = 0. (16) 

In the present study, no nodal lines are found in < 
J < IOTl. Nodal lines are found in 12?i < J < 20fi. and 
form a simple topology of closed loops (rings) without 
knots. Their geometrical distributions seem to become 
more complicated as the value of J becomes larger. It is 
seen in FigEJthat some of the nodal lines are missing in 
the low resolution calculation but this result is due to the 
simple linear approximation, Ea. (|15(l . for the complicated 
geometrical distribution of the nodal lines. 

Symmetries seen in Figs El and are expressed by 

(R(a,f3,"/)) = (R(j,P,a)}* (17) 
= (R(tt - a, 7T - p, j))*=(R(a, 7T - A 7T - 7 ))* (18) 
= (R(a + 2ir,/3,a)) = (£(«,/?, 7 + 2vr)}. (19) 

and are attributed to the properties possessed by the 
principal-axis cranked HFB states ^7| such as the sig- 
nature. As a result of these symmetries (especially, Eqs. 
(fTTjl and (CHJ), there are four identical sets of nodal lines 
in Figs El and |3 In other words, there is essentially only 
one nodal line (ring) out of the four shown in FigsEJand 
EfA), while there are sixteen nodal rings in FigsEjB) but 
just four of them are essential. Two out of these four es- 
sential rings form a concentric circles, as shown in the 
second portion of Fig^Jb). (Note the symmetry on the 
boundary is given in Eqs. 119(1 .) Although these rings look 
to be completely two-dimensional, that is, confined in the 
constant (3 plane, one can tell their three-dimensional 
structure with a careful look at the figures. 



Let us now explain the details of the new approach to 
determine the sign a{Vt) below. First of all, the third 
order derivatives of the norm overlap kernels are chosen 
to be a measure of the "smoothness" of the kernels. In 
fact, such a quantity is practically easily computed in 
our calculations where the norm overlap kernel and its 
first order derivatives are already obtained. Using the 
Taylor expansions for N(fl + Afi) and dN ^Q An ' > , the 
third derivative between two points f2 and fl' — f2+ A51 is 
approximately obtained up to an order of 0(Af2), which 
is denoted here as D 3 (Q',Q), 

6 i dN dN 1 

D ^ = (A^W^+a^} < 20 > 

The above expression corresponds to a case when two 
points arc displaced along the a coordinate, i.e., = 
(a,/3, 7) and fi' = (a + Aa,/?, 7). Similar expressions 
hold for a pair of adjacent points displaced in the (3 or 7 
direction. Note that D3 is a complex quantity in contrast 
to the residue 89, which is a real quantity. An advantage 
to extend the residue 89 to the complex D3 is that the 
continuity is checked not only by 0(0) but also by r(O). 
-D3 is particularly useful in the vicinity of the nodal lines 
where lnr(£!) diverges so quickly. In this case, 89 can 
also change rapidly to return the large value as if the 
boundary of domains appeared to be detected. With 
all the available information on the continuity both in 
the norm and argument, the boundary between domains 
D n and D n+ \ is determined with more certainty even 
in the vicinity of the nodal lines. Next, starting from 

= (0,0,0), where a = 1 is assumed, we determine one 
by one the sign of a for the adjacent mesh points so that 
I-D3I takes a smaller value. As we proceed, conflicting 
assignments in sign are seen for a point fl = (a,/3, 7). 
In other words, when D3 is evaluated along two different 
links, say links between fl and + Ail = (a ± Aa, /?, 7), 
the signs at 17 are sometimes inconsistent. In this case, 
we use "reliability" to decide which sign assignment has 
priority. Although there is no unique way to define the 
"reliability", our choice for the definition is the ratio of 
I-D3I between the sign-flipped and unflipped cases (the 
smaller value should be in the denominator): larger val- 
ues of this ratio indicate more reliability. This proce- 
dure is performed for all the links and optimisation on 
£>3 is carried out. The whole optimisation procedure is 
repeated many times until the average reliability is con- 
verged. The essential point in this new method is that 
we attempt to proceed with the analytical continuation 
through the only links where a reliable determination of 
the relative sign seems possible. With this requirement, 
many problematic points like those near the nodal lines 
are circumvented or isolated from the analytical contin- 
uation procedure. 

This new method is quite powerful especially when 
one attempts to determine the sign of norm overlap ker- 
nels in the presence of the nodal lines. An advantage 
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FIG. 5: Probability distributions for J = 12h - 20ft with the 
modified phase determination method. The upper (lower) 
panel shows the profile for even-/ (odd-/) components. 



of the method is that we do not need to find the nodal 
lines explicitly. To confirm the advantage, we analysed 
where the disagreements in the sign determination hap- 
pen between the previous naive method (in which the 
signs are determined one by one in a scheduled simple 
order along straight lines) and the present new method. 
We then found that the discrepancies are concentrated in 
the vicinity of regions where nodal lines exist. The im- 
provement in the sign determination by the new method 
is clearly demonstrated in FigQ as the recovery of the 
positive definitcness of the probability in the odd-/ distri- 
butions of W 1 . A disadvantage of this new method is that 
it is more time-consuming than the previous method. 
However, we found that the computation time has been 
greatly improved with the combination of the two meth- 



ods: the initial sign determination is made by the first 
naive approach and subsequently the new optimisation 
method is carried out. 

In the case of 170 Dy, we have checked that this newly 
developed method improves the sign determination up 
to J = 14ft. Although several doubtful sign assignments 
start to occur at J = 16h even for the higher resolution, 
these potential mistakes do not seem to affect positive 
definiteness in probabilities (|C^/| 2 > 0) and the sum 
rule, 



W 1 = 1, 



(21) 



up to J = 20Ti. In Fig|5J the results of the modified 
probability distributions are plotted. The integrals of 
each probability distribution curve, i.e., the probability 
sum rule, are presented in Table HJ W 1 satisfies the sum 
rule quite well up to / = 34?i. The problem of nega- 
tive probability does not appear until / = 3lh. Beyond 
/ = 3lh for the odd-/ components (and beyond / = 34?i 
for the even-/), all W 1 become negative. The reason for 
this behaviour is attributed to the discretisation approx- 
imation in the integration. With Aa = A7 = 5° and 
A/3 = 1°, the integration in angular momentum projec- 
tion manages to be valid up to 



2h 
Aa 



2h 2h 
A/3 + A7 



32ft. 



(22) 



The divergence in the curves of J = 12ft (beyond / > 34ft 
for the even-/ graph and / > 33ft for odd-/) comes from 
the same reason. To overcome even these difficulties, 
a finer mesh size should be employed. Such a line of 
investigation is now in progress. 

In conclusion, we have studied the distribution of 
zeros of norm overlap kernels of cranked HFB states. 
Such zeros form one dimensional structure in the three- 
dimensional space defined by the Euler angles, so that 
we call it a "nodal line" . The existence of nodal lines 
is numerically demonstrated for the first time, and we 
have found they emerge as spin increases and their struc- 
ture becomes more complicated as angular momentum 



J = (Ji)hpb 


12h 


Uh 


16h 


18h 


20h 


p 


0.998 


0.991 


0.981 


0.965 


0.939 



TABLE I: The total probability sum P = ^ W shown in FigEJ 

K35R 



becomes larger (J > lOTi). We have shown that it is 
important to know the distributions in order to correct 
wrong assignments for the sign of norm overlap kernels. 
The correction is essential for angular momentum pro- 
jection to be performed with high accuracy for cranked 
mean-field states. 
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